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Abstract 

Within the framework of the extended teleparallel gravity, a new class of boson 
stars has recently been constructed by introducing the nonminimal coupling of the 
scalar field to the torsion scalar. An interesting feature of these static, spherical, 
self-gravitating configurations of the massive complex scalar field is their central 
region with outwardly increasing energy density, surrounded by a thick shell within 
which the joining with the usual asymptotically Schwarzschild tail takes place. In 
this work we extend the original model with the U( 1) gauge field and we find 
that the combined effect of the charge and coupling of the field to torsion leads 
to a significant increase of the maximal mass and the particle number that can be 
supported against gravity. We also show that some charged configurations preserve 
the property of having the outwardly increasing energy density over the central 
region, regardless of the fact that charging the configurations affects the anisotropy 
of the pressures in the opposite way relative to that of the field-to-torsion coupling 
terms. 


1 Introduction 

Our motivation for considering the nonminimal coupling of the scalar field to the torsion 
scalar, and boson stars in particular, comes from the recently revived interest in the 
theories of gravity based on torsion, rather than curvature in a- It has been known for 
a long time that formulating the gravitational action with the torsion scalar T, instead 
of the curvature scalar R as in the standard theory, leads to the theory of gravity that 
is equivalent to general relativity (GR) and that is known as the teleparallel equivalent 
of GR (TEGR). This theory uses the curvature-free Weitzenbock connection and its 
dynamical variable is the tetrad—a set of four orthonormal fields that select the local 
Lorentz frame at every spacetime point. The extensions of the TEGR are particularly 

* Email: sasa.ilijic@fer.hr 


1 



appealing because they are in some aspects simpler than the similar extensions of GR. 
For example, the f(T) theory, which could be compared with the f(R) extension of 
GR, leads to equations of motion that are of the second order, whereas f(R) leads to 
equations of the fourth order. This is due to the fact that the torsion scalar involves 
only the first order derivatives of the tetrad, while the curvature scalar involves the 
second order derivatives of the spacetime metric. However, there are also complications 
associated with the extensions of the TEGR that are not present in the similar extensions 
of the standard curvature-based theory. Within the TEGR there is the six-fold infinity of 
acceptable tetrads, related among themselves by rotations and by Lorentz boosts, which 
reflects the Lorentz invariance of the theory. In the simplest extensions of the theory, 
such as the f(T) theory [31 0], or theories involving nonminimal coupling of the scalar 
field to the torsion scalar T EM, the Lorentz invariance of the resulting equations of 
motion is being lost. The equations of motion must be engaged to single out the right 
tetrad, i.e. one must deal with the six degrees of freedom of the tetrad (related to local 
rotations and Lorentz boosts), for which there is no universal procedure. Regardless of 
these difficulties, in several special cases that include the static spherical symmetry and 
homogeneous isotropic geometry, the preferred tetrads leading to physically meaningful 
equations of motion have been found [8j and the solutions that are distinct from their 
extended GR counterparts could be constructed and studied. 

In the static spherical symmetry, boson stars are a good starting point to study 
the effects of newly introduced terms, such as those due to scalar-to-torsion coupling, 
because they are arguably the simplest mathematical models of static self-gravitating 
objects with the underlying matter model derived directly from the field theory. Within 
the standard GR, boson stars were first constructed by Kaup [9] and by Ruffini and 
Bonazzola ra as self-gravitating configurations of the massive complex scalar field. The 
basic model was soon extended with the inclusion of the field self interaction HU , electrical 
charge ra, nonminimal coupling |13j . etc. This broad subject has been reviewed in 
ns nans nzi . Boson stars with nonminimal coupling of the scalar field to the curvature 
scalar revealed exotic properties such as configurations with negative principal pressures 
in the interior of the star [18] or the formation of the photon spheres in the strong 
gravitational lensing regime [19]. Nonminimal coupling of the scalar field to the torsion 
scalar revealed another interesting property—the outwardly increasing energy density 
over the interior of the star [20]. Such behavior has not been encountered in any other 
stellar model with energy-momentum derived from a simple field theory. In this work we 
intend to study this feature more closely and also inspect the behavior of the anisotropy 
of the principal pressures, which is a mechanism behind many of the exotic spherically 
symmetric solutions in general relativity, e.g. wormholes [21] or gravastars (stars with 
negative pressures in the center) [22]. We also extend the torsion-coupled model by 
adding the local 17(1) gauge invariance, or the electric charge, a component that has not 
yet been considered in this context. As the static scalar field and the electrical field are 
both sources of the pressure anisotropy in the static spherically symmetric setting, we 
are interested to investigate how the outwardly increasing energy density obtained by 
scalar-to-torsion coupling will respond to the inclusion of the electrical charge. 

The paper is organized as follows. A necessary minimum of the technical background 
on the teleparallel gravity, mainly restricted to establishing the notation, is given in Sec. [2] 
Static spherical symmetry is introduced in Sec. [31 and in Sec. Q] we construct a number 
of specific boson star configurations and examine the behavior of the energy-momentum 
components. We discuss our results in Sec.[5j Some additional technical details are given 
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in the appendix. Geometrized units (c = 1 = Gn) are used in the paper. 


2 Charged scalar field coupled to torsion 


The tetrad fields are denoted by h a M , where the latin index runs over the Lorentz frame 
coordinates, and the greek index runs over the usual spacetime coordinates. The tetrad 
satisfies the orthonormality relations g ab = h a fJ 'hb u g fJ , u , g^u = h a u r} abr h a and 
ff jij = 6&, where g a b = diag(—, +, +, +) is the metric in the Lorentz frame, and g /w 
is the spacetime metric tensor. Torsion-based teleparallel gravity uses the curvatureless 
Weitzenbock connection that is defined in terms of h a [L as 

r a 01 = h a a d,h a 0 . (i) 

The torsion tensor is defined as 

T ex _ pa pa /o\ 

M = 1 7/5 _ 1 \ Z ) 

and the torsion scalar as 

f = \f aM f a ^ + l -f a01 f^ a - t Q7 q t%. (3) 

The essential property of T is its relation to the curvature scalar R , 

R = -f~j^(hf^ u ), (4) 

where h = det [h a v \ = \J— det[g a0 ], showing that the two scalars differ only by a total 
divergence term. This implies that while the action of the TEGR, 


Stegr — 


dx 4 h 



(5) 


is not invariant with respect to Lorentz transformations of the tetrad (i.e. the choice of 
the tetrad), the resulting equations of motion are equivalent to those of GR because the 
variation of the action is not affected by the divergence term. 

Here, we consider the action S — f C h d A x with the Lagrangian density 

£ = (-T) - + D^DAT) - » 2 <t>*<t>, ( 6 ) 

which involves the massive complex scalar field (j) nonminimally coupled to T and the 
gauge field A^. The quantities k = 87t, £ and e are the coupling constants and fi is the 
field mass. = d fJj A u — d u A fl is the usual field-strength tensor and D /t = d tl — ieH /t is the 
gauge-covariant derivative. Variation of the action with respect to the tetrad leads to the 
equation of motion that can be written in the form of the Einstein equation, G^ = hT /Jiy , 
where G^ = R IIU — | Rg /W is the Einstein tensor and is the energy-momentum tensor. 
The energy-momentum tensor can be written as 

rffl - 2^ u a d a (r^) 

^ 1 + 2k£(f)*(f) 
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where 


T$ = F^F u a + (D „(/))* D„(j) + D^D^)* 

- <V (jFafiF'** + + D a (t>(D p <t>)*) + m 2 0*<^ (8) 

is the familiar energy-momentum tensor that one obtains in the case of the minimal 
coupling of the scalar field, either to R in the standard theory, or to T within the torsion- 
based TEGR. The denominator of (J7J) is the same as in the case of nonminimal coupling 
of the scalar held to the curvature scalar, while the second term in the numerator reveals 
the difference. The tensor S a ^ f that appears in this term is known as the modified torsion 
tensor and it can be written in terms of the torsion tensor (J2]) as 

= ^ {Tpar/ T -^7/3a T -^"a/37) T 9 a /3 Ty 9 a 7 F3 • ( 9 ) 

The variation of the action with respect to the scalar held gives 

—dJhD^cj)) = ieAnD^cj) + (FT + m 2 + 2A 0*0)0, (10) 

h 

and the variation with respect to the gauge held gives 

\d,(hFn = ie((ZT»0* - (D^Ycf). (11) 

The Lagrangin density dHJ) also implies the conserved current = i((_D /i 0)*0 — _D M 00*) 
with the particle number N and the charge Q = eN as the conserved quantities. 


3 Spherical symmetry and anisotropy of pressures 

In order to study the static spherically symmetric configurations of the charged scalar 
held we adopt the usual coordinates (t, r, ’d,(p), where r is the area radius, and we write 
the line element as 


ds 2 = _ e 2 *h)/ 2 ( r ) dt 2 + /~ 2 (r) dr 2 + r 2 dfl 2 , (12) 

where dfl 2 is the line element on the unit 2-sphere, and $ and / are the r-dependent 
metric prohle functions. For the scalar held we take the time-stationary Ansatz 

0 = 0( r ) e --\ (13) 

and from here on 0 will denote only the real r-dependent part, and 0' will denote its 
radial derivative. For the gauge held we assume 

A* = <W r )> ( 14 ) 

i.e. the held has only the radially dependent electrostatic potential -0. In the static 
spherical symmetry the energy-momentum tensor has the form T£ = diag(— p,p, q, q). 
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where p is the energy density and p and q are the radial and the transverse pressures. In 
the case of minimal coupling (£ = 0) these quantities follow from (JH]) as 

P (0) = + /V 2 + r 2 e- 2 > + e-0) 2 + (15) 

P (0) = -h 2 0 2 + f 2( P' 2 + f~ 2 e~ 2 *{u + e-0) 2 - ^e“ 2< V 2 , (16) 

g (°) = _ y20/2 + f-2 e -2*( u + g 0)2 + (17) 

In order to include the contribution to the energy-momentum tensor due to nonminimal 
coupling to the torsion scalar one needs the modified torsion tensor. The tetrad must be 
chosen first and we adopt the form 

0 0 0 \ 

Z” 1 sill'd cos <p r sin'dcos cp — r ski'd sin tp , . 

/ _1 sin d sin (p r sin r) cos cp r sin 'd cos ip ’ ^ ' 

/ _1 cos'd —r sind 0 J 

which has proven to lead to physically meaningful equations of motion in the context of 
static spherical symmetry [23]. A lengthy but straightforward calculation, which involves 
the computation of the Weitzenbock connection (jT]) , the torsion tensor (|2j), and the mod¬ 
ified torsion tensor (|9jl, leads to the following expressions for the energy density and the 
pressures: 


/e*/ 

0 

0 

V 0 


p (0) + 8 UU_ - 1 )(M>'/r 

1 + 2 k£<f> 2 

p( o) 

1 + 2k Z0 2 ’ 

q(o) _ 4£f(f — 1 + rf + rf^')(j)f/r 

1 + 2k£(j) 2 


(19) 

( 20 ) 
( 21 ) 


where p^°\ p^ and q^ are given by fl 151) - fl 1 7D . The anisotropy of pressures, which we 
dehne as the difference between the transverse and the radial pressure, follows as 


-2/V 2 + e- 2, V 2 - 4 £f(f - 1 + rf + rf&^f/r 
1 + 2 k£cj) 2 


( 22 ) 


It is easy to see that in the absence of the electric charge (e = 0, ip — 0) and minimal 
coupling (£ = 0 ), the anisotropy reduces to q — p = —2 f 2 f 2 which is negative, except at 
r = 0 and as r —y oo where it is expected to vanish. 

It is also straightforward to express the equations of motion (TT01) and (fTTl) in terms of 
our metric and held Ansatze, but as the resulting expressions are somewhat lengthy, we 
only show the expression used to compute the particle number, 

/ poo 

d 3 x = 87 t J r 2 e _$ / _ 2 0 2 (a; + eip) dr (23) 


(Q = eN is the charge). The total mass of the configuration is given by the familiar 
expression M = An J 0 °° r 2 p(r) dr, where p is the energy density given by ffl9]l . 
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Figure 1: Mass (solid lines in units MpJn) and particle number (dashed lines in units 
Mpj/// 2 ) of boson stars with parameters given in Table |T] and for a range of values of 
(Jo = /c _1 / 2 0(0). Critical configurations are indicated with vertical lines and letters. 


4 Properties of the charged torsion-boson stars 

As the solutions to the equations of motion are not accessible in a closed form, we proceed 
with constructing them numerically. Solutions can be obtained by solving a boundary 
value problem for the metric profile functions $ and /, the fields 0 and 0, and the 
eigenvalue u>, with the boundary conditions specified at the center of symmetry (r = 0) 
and at the spatial infinity (r = oo). At r = 0 we require 0 = 0 O , < h / = 0 and / = 1, 
while at r = oo we require 0 = 0 = 0. These boundary conditions, together with 
the Einstein equation which provides three independent differential equations, and the 
equations of motion for the fields (flOl) and (fTTf . fully determine the solution. The value 
of 0o (scalar held amplitude at r = 0) that enters through the boundary condition is 
used to parametrize the family of solutions corresponding to the fixed values of the held 
mass fi and the coupling constants e (electric charge) and £ (nonminimal coupling of the 
scalar held to the torsion scalar). For a chosen set of the values of ^ e and £ we can 
compute a spectrum of solutions starting from a very small value of 0o- (We have used 
/r = 1 with no loss of generality, see [12] for the rescaling scheme that can be applied. 
We have also limited our analysis to solutions without nodes in the scalar held.) As 0o 
increases, the particle number and the mass of the configuration increase as well, until they 
simultaneously reach the respective maxima. The configurations with maximal N and M 
are known as critical because it in the case of the minimal coupling of the held they mark 
the transition from the dynamically stable to the dynamically unstable configurations. 
This behavior of N and M is shown in Fig. [1] for four sets of parameters. The critical 
configurations are indicated with letters A to D, and their most important parameters 
are listed in Table [EJ As it is clear from the data in the table, we have constructed the 
configurations corresponding to the minimally coupled electrically neutral held (model A), 
minimally coupled electrically charged held (model B), nonminimally coupled electrically 
neutral held (model C), and the critical configuration with the nonminimally coupled 
electrically charged held (model D). One can see that introducing the charge alone makes 
it possible to support a larger number of particles against gravity and obtain a stable 
boson star of larger mass m- The same holds for introducing of the held-to-torsion 
coupling, while the combined effect of charge and torsion coupling leads to a still larger 
increase in the critical N and M. 
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Table 1: Parameters of the four critical configurations of the boson stars. Models C and 
D include nomninimal coupling of the scalar field to the torsion scalar. 


Model 

h 

e 

£ 

v 0 

CO 

N/iM^/p 2 ) 

M/(Ml/n) 

A 

1 

0 

0 

0.271 

0.853 

0.653 

0.633 

B 

1 

3 

0 

0.184 

0.947 

1.330 

1.316 

C 

1 

0 

8 

0.575 

0.685 

1.008 

0.927 

D 

1 

3 

8 

0.343 

0.890 

2.184 

2.119 


The model A represents the critical configuration of the minimally coupled electrically 
neutral field. Since the earliest works on the boson stars, this has been the common 
reference configuration when considering extensions of the simple model. The energy 
density and the pressures are shown in the upper left panel of Fig. [2j As is expected on 
the basis of (1221) . the anisotropy of the pressures is negative everywhere except at r = 0 
where it must vanish by virtue of symmetry and at spatial infinity where both p and q 
vanish asymptotically. The transverse pressure crosses zero at finite r and asymptotically 
approaches zero from the negative regime as r —> oo. 

The model B includes the electric charge. The value of the coupling constant e = 3 
is chosen close to the maximal value that the configuration can support. The upper 
right panel of Fig. [2] shows that, on the scale set by the central energy density of the 
configuration, the anisotropy of pressures is of small amplitude. However, this in no way 
means that its role is unimportant. The anisotropy of pressures in the interior is negative, 
which is the characteristic feature of the minimally coupled field, while in the asymptotic 
tail of the configuration it becomes positive. This reflects the fact that, as r —)■ oo, the 
metric is asymptotically Reissner-Nordstrom, which implies p = — p = q = Q 2 /(87ir 4 ), Q 
being the electric charge. In this model the radial pressure changes sign and is negative 
in the asymptotic tail, which is exactly the opposite of what we had in model A, and is 
a characteristic of electrically charged configurations. The coupling constant e (charge) 
can be further increased, but as one approaches the maximal amount of charge the boson 
star can support, the anisotropy becomes positive over the whole configuration. 

Model C includes the nonminimal coupling of the scalar field to the torsion scalar, 
but no electric charge. The value of the field-to-torsion coupling constant £ = 8 has been 
chosen large enough for the configuration to develop the outwardly increasing energy 
density over the interior. The behavior of pressure anisotropy is exactly the opposite 
of what we had in model B. In the interior we have the positive anisotropy and in the 
asymptotic tail it becomes negative. The transverse pressure has a zero crossing, while 
the radial pressure remains positive. If one used a larger value of £ the above properties 
would only be further amplified, but would not change qualitatively. 

In model D we include both the electric charge and the nonminimal coupling of the 
scalar field to the torsion scalar. The values of the coupling constants are those that 
were used individually in models B and C. The combined effect of charge and torsion 
preserves the outwardly increasing energy density over the interior and results in the 
positive anisotropy of pressures throughout the boson star. The radial pressure has the 
zero crossing. This is exactly the opposite of what one finds in a simple boson star. As 
in model B, the coupling constant e could be further increased, but as it approaches 
the maximal (critical) value of the charge that can be supported by the boson star, the 
outwardly increasing energy density changes to outwardly decreasing, which means that 
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Figure 2: Energy-momentum components of the critical configurations A-D (see Table Q] 
for parameters): energy density (solid thin line) relative to the central energy density po , 
radial and transverse pressure (solid and dashed pair of thin lines) relative to po , and the 
anisotropy (gray thick line) relative to p 0 and scaled up by a factor for clarity. 
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the configuration becomes qualitatively similar to model B. 


5 Discussion and conclusions 

In the static spherically symmetric setting, the scalar field produces the negative anisotropy 
of principal pressures (defined as the transverse minus the radial pressure), while the ra¬ 
dial electrostatic field gives a positive contribution to the pressure anisotropy. These two 
components contribute in the opposite directions, so one can speak of the scalar and 
of the electrostatic sign of the anisotropy. The matter model that we considered also 
includes the coupling of the scalar field to the torsion scalar, which is the third mech¬ 
anism contributing to the pressure anisotropy, and also the component giving raise to 
the outwardly increasing energy density in some configurations. From the expression for 
the anisotropy (1221) it is not obvious in what way the torsion-coupling terms affect the 
anisotropy. In addition, there is no direct way to isolate the mechanism responsible of 
generating the outwardly increasing energy density. We therefore obtained numerical 
evidence that shows the trends that take place when the three sources of anisotropy are 
active simultaneously. A small amount of the charge added to a simple boson star changes 
the sign of the anisotropy from the negative (scalar) to positive (electromagnetic) regime 


























only in the outer layer, i.e. in the asymptotic tail, while as the amount of charge ap¬ 
proaches the critical value beyond which the charged star does not exist, the anisotropy 
becomes positive throughout the configuration. The nonminimal coupling of the field to 
torsion, but without charge, changes the sign of the anisotropy into positive (electromag¬ 
netic) only in the interior part of the configuration, while the anisotropy remains negative 
(scalar) in the asymptotic tail. If both charge and torsion coupling are present simulta¬ 
neously, provided that field-to-torsion coupling and charge are not too weak, we find that 
the anisotropy is positive throughout the configuration and also the outwardly increasing 
energy density over the interior of the boson star. The energy density returns to the usual 
outwardly decreasing behavior only as the charge approaches the maximal amount that 
can be supported. The unusual property of outwardly increasing energy density is here 
shown to be present in a wider range of physical situations than it was known previously 
m Therefore the torsion-coupled scalar field confirms itself as an interesting candidate 
for the modeling self-gravitating structures. It is also worth noting that torsion coupling 
of the scalar field helped support a larger number of particles against gravity, which is a 
trend similar to that of the charge [12] or of the quartic field self interaction [Ill- 

Before closing, we should not forget to mention that the matter model we considered 
was constructed within the framework of the extended teleparallel gravity, which means 
that, in order to work with physically meaningful equations of motion, a tetrad suitable 
for the assumed geometry had to be chosen first. As we considered the static spherical 
symmetry, we could use the good tetrad provided by [Hj. However, applications of this 
matter model outside of the static spherical symmetry or homogeneous isotropic geometry 
could be hampered by the difficulties in finding the suitable tetrad. For example, the 
tetrad compatible with the time-dependent spherically symmetric geometry seems to be 
still out of reach. 

Acknowledgments: This work is supported by the University of Zagreb Grant No. VIF2014- 
PP2-11. Partial support comes from “NewCompStar,” COST Grant No. MP1304. 

A Explicit equations of motion and the numerical 
technique 

For completeness here we give the explicit form of some expressions that are, in order to 
keep the presentation compact, not included in the main body of the paper. 

Adopting the metric (fT2jh the Ricci or curvature scalar reads 

R = ^ - 2 f 2 - y (/'(4 + 3r$') + rf ") - ^((1 - r<F') 2 + r 2 $"), (24) 

and the three nontrivial components of the Einstein tensor G^ u = R^ u — \R8f; are 

G f t = ^(-l + / 2 + 2r/f), (25) 

G\ = - 1 + / 2 (1 + 2r$') + 2 rff), (26) 

G 9 # = G% = f 2 + - (/'(2 + 3r&) + rf") + £ ($' + r$' 2 + r<f>"). (27) 

The Einstein equations, G^ u = are obtained by combining (1251) - (1271), and the 

components of the energy-momentum tensor, = diag(— p,p, g, q), that are given by 
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m-m an d (dSD — (gZD - Adopting the tetrad (fTSlh the torsion scalar (J3]) , which is related 
to the curvature scalar by (j4]), reads 

+ /«')). ( 28 ) 

The equation of motion for the scalar held (fTUj) can be written as 


e 2 ®(p(u} + e'l/i) 2 4>'{2rf + /(2 + r<h')) (pp 2 

¥ rf + 

2^(1 - f)(l ~ f ~ 2r(f + /$')) 


where the terms that are specific to the scalar-to-torsion coupling are grouped in the 
lower line. The only nontrivial component of the equation of motion for the gauge held 
(ITT]) is its time component. It does not involve terms due to scalar-to-torsion coupling 
and it reads 

;// 2e4> 2 {uj + e-0) -0'(2 — rd*') 

I 2 r 

Our procedure for solving the system of nonlinear ordinary differential equations for the 
four unknown functions, /, <L, 0 an d 0, subject to the appropriate boundary conditions 
(see Sec. ED, consists of changing the radial variable r —> x = r/(r + 1) in order to 
formulate the problem on the compact domain and of solving the resulting boundary 
value problem (BVP) numerically. The frequency u that appears in the harmonic Ansatz 
(USD f° r the scalar held has the role of the eigenvalue and is obtained for every converged 
solution. We used the collocation algorithm implemented in the numerical code COLSYS 
|24j . In order to double-check the validity of the numerical solutions, all of them were 
confirmed by the independent initial value integration that can be easily carried out once 
the eigenvalue c o has been determined by solving the BVP. 
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